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ABSTRACT 

Optical resonances in ID photonic crystal microcavities are investigated numerically using finite-element light 
scattering and eigenmode solvers. The results are validated by comparison to experimental and theoretical 
findings from the literature. The influence of nanometer-scale geometry variations on the resonator performance 
is studied. Limiting factors to ultra-high Q-factor performance are identified. 
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1. INTRODUCTION 

Integration of optical components with electronic circuits promises to provide increased functionality. Building 
blocks of such integrated optical devices are waveguides, modulators, switches, and others. Due to their high 
quality factor (Q-factor) and small modal volume, one dimensional photonic crystal (PhC) microcavities in 
silicon-on-insulator are candidates for integrated all-optical switching devices. High Q-factors have been reported 
for such devices by Zain et al} Further realized designs of optical microcavities include 2D photonic crystal 
resonators--- suspended air-bridge cavitiepl and circular grating resonators.!^ 

In this contribution we numerically investigate high-Q resonances in PhC microcavities, where we partially 
review results of an earlier contribution.'^ We analyse the photonic devices through numerical simulations of 
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Figure 1. Geometry of a PhC microcavity, discrctizcd with a prismatoidal mcsli. Dark grey: Silica substrate, light grey: 
Silicon waveguide, perforated by air cones (side- wall angle of 85 degrees). The structure is surrounded by air, substrate 
and waveguides extend to infinity. 



Figure 2. Electric field intensity distribution in the microcavity close to a resonance. Cross-section at z = h/2 above the 
substrate through the 3D solution in pseudocolor representation (blue/dark: low intensity, red/bright: high intensity). 



light propagation using a time-harmonic finite-element method (FEM) light scattering solver. The method is 
validated by comparison to experimental and theoretical results, and the effect of small geometry variations 
on device performance is investigated. We further use an eigenmode solver for direct computation of resonance 
wavelengths and Q-factors. The light scattering solver, the eigenmode solver as well as a propagation mode solver 
used for computing incident waveguide modes are parts of the commercial FEM programme package JCMsuite, 
which is developed by JCMwave and ZIB.^ The solvers enable fast and accurate simulations of main properties 
of optical microcavities. 

2. INVESTIGATED SETUP 

The investigated setup consists of a silicon waveguide with square shape which is perforated with a number of 
cylindrical air pores (PhC) and is supported by a silica substrate. The cavity essentially is formed by a missing 
central pore. Diameters and positions of the air pores are fine-tuned for high quality and transmission.^ For the 
FEM simulations we discretize this geometry by prismatoidal elements. Figure [T] shows a part of the spatial mesh. 
A typical computational domain has a size of 1.6 /im x 6 /imx 1.6 /im corresponding to roughly 3.5 x 13 x 3.5 ~ 160 
cubic wavelengths (in Si). The model extends waveguide, substrate and surrounding air to infinity by applying 
transparent boundary conditions to the corresponding faces. The material relative permittivities are = 3.48^ 
(Si) and — 1.44^ (Si02), the waveguide width and height are w = 520 nm and h — 340 nm. Air holes 
with diameters di = 130 nm, ^2 = 160 nm, — 185 nm, ^4 = 181 nm, d^ = 181 nm (holes with d^ repeat 
= 2 . . . 4 times) , are positioned symmetrically around the cavity center with center to cavity-center distances 
±Pi: pi — L/2 + di/2 (where L is the cavity length, L = 422 . . .425nm), p2 — pi + 300 nm, ^3=^2+ 315 nm, 
P4 = P3 + 325 nm, p^ = P4 + 352 nm, pn+i = Pn + 370 nm. This setup has been investigated experimentally and 
theoretically by Velha et al for L = 425 nm.'^ We have extended the investigation to a fine sampling of cavity 
lengths L on a nanometer scale, and we investigate pore sidewall angles different from 90 degrees. 

For simulating transmission through the microcavity at specific wavelength Aq, we first compute the funda- 
mental propagation mode of the waveguide at Aq using a FEM propagating mode solver.lSl The obtained mode 
field is applied as input data to one of the boundaries of the 3D computational domain (left boundary in Fig. [1} , 
such that the mode propagates in direction of the center of the waveguide. We then compute the scattered 
light field in the setup corresponding to this excitation using higher-order finite-elements. In post-processes we 
extract energy fluxes through interfaces and field distributions in several cross-sections from the 3D light field 
distribution. Figure [5] visualizes the computed electric field intensity at an incident wavelength close to the 
cavity's resonance wavelength. 

3. NUMERICAL RESULTS 

3.1 Validation of FEM results 

In order to compare results obtained with our method to results obtained with a dedicated Fourier modal method 
(RCWAj^'^ and to experimental result^ we have performed wavelength scans over cavity resonances for cavity 
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Figure 3. Transmission spectra for PhC cavities with L — 425 nm and A'^ — 2/3/4 on a logarithmic scale. 

lengths L — 425 nm.- Figure [3] shows the corresponding transmission spectra for cavities with a total nmxiber 
of 12/14/16 PhC holes, resp. for TV = 2/3/4. As expected, for all three cases we find clear resonances in the 
transmission spectra. All resonances are centered around a peak wavelength of about Xq^fem = 1593 nm. From 
the ratio of Aq to the widths of Lorentz fits to the distributions we find Q-factors of Qfem — 880/2630/7580 
for N — 2/3/4. The simulated peak wavelengths correspond very well to experimental results of Velha et al 
(Ao,exp ~ 1599/1592. 7/1592. 7nm for TV = 2/3/4) and to numerical results of Velha et al {Xo,rcwa « 1600 nm). 
The Q-factors obtained with our method are slightly below the experimental and numerical values from Velha 
et al {Qexp = 1100/2700/8900, Qrcwa = 1270/3200/9600). 

To further validate our results we have performed a convergence analysis showing that with moderate com- 
putational effort in the range of minutes of computation time for a single frequency we reach accuracies of the 
resonance wavelength of better than 0.1%, and of high Q-factors^ in the few percent range. 

3.2 Investigation of geometry variations on a nanometer- scale 

We have performed numerical experiments with slight variations of some geometrical parameters.- This allows 
to quantify the influence of fabrication tolerances on microcavity properties. Due to spatial adaptivity of the 
used unstructured meshes, FEM allows to investigate influences of very small variations of the geometry at high 
accuracy and without additional computational costs. 

In a first set of simulations we have performed wavelength scans for PhC cavities with 12 air pores (N — 2) 
for cavity lengths of L = 422, 423, 424, 425 nm. Figure U (left) shows the corresponding transmission spectra. As 
can be seen from these spectra and from the fitted distributions, in the investigated regime, a shift of L shifts 
the resonance wavelength Aq approximately linearly with 6X0/6L w 0.8. Also, a decrease of L by about one 
percent increases Q by about 4% in the investigated regime. This shows that fabrication tolerances in the range 
of few nanometers can well explain significant deviations of device performance. The slight shift of resonance 
wavelength Aq with respect to the data displayed in Fig. [3] of about 0.02% is due the lower numerical resolution 
of the simulations displayed in Fig. 2] (left). 

In a second set of experiments we have varied the sidewall angle (SWA) of all air pores between 84 and 90 
degrees (90 deg. corresponds to perpendicular SWA, resp. cylindrical shape). For a better comparability of the 
results, for SWA < 90 degrees the upper radius of the air pore was increased, the lower radius was decreased, and 
the radius at central waveguide height was kept constant. Figure 3] (right) shows corresponding results: While 
the influence of pore SWA on resonance wavelength can be neglected in the investigated regime, the influence of 
pore SWA on Q-factor is signiflcant: A deviation of 5 degrees from vertical sidewalls decreases the Q-factor by 
about 10 percent. 



3.3 Direct simulation of cavity resonances 

Alternatively to simulating the light scattering response of the device, we can also use an eigenmode solver for 
directly simulating the resonance properties (resonance wavelength and Q- factor): Given the geometrical setup 
as described in Section [2l one computes an electric field distribution E and a complex eigenfrequency uj which 
satisfy Maxwell's time-harmonic wave equation 

V X ^"^V X E = uj^eE 

without sources; electric permittivity and magnetic permeability are denoted by s and fj,, respectively. Transpar- 
ent boundary conditions take into account the specific geometry of the problem where waveguides are modelled 



□ L = 425 nm 
fit, Q = 7537 




Wavelength [nm] Sidewall angle [deg] 

Figure 4. Left: Transmission spectra for PhC cavities with L = 422... 425 nm and TV = 2. Right: Dependence of 
relative changes of Q-factor and resonance wavelength with sidewall angle (for PhC cavities with L = 425 nm and N — 2). 
Circles correspond to relative change of the Q-factor, SQ — (Qsqa — Q9o)/Q90- Squares correspond to relative change 
of resonance wavelength, SX = {Xo,swa — Ao,9o)/Ao,9o, magnified by a factor of 1000. Q90 and Ao,9o are the values for 
SWA=90deg. 
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Figure 5. Convergence of the Q-factor (left) and of the resonance wavelength (right) obtained with a FEM eigenmode 
solver for polynomial finite element degree p = 2,3, 4. 
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Figure 6. Computational costs corresponding to the simulation results displayed in Figure[5] Left; Number of unknowns 
in the corresponding system of linear equations, A''. Right: Computation time on a workstation with four double-core 
AMD Opteron CPU's at 2.8 GHz has been used. Maximum memory consumption was about 32 GB of RAM (for N=4, 
p = 4). 



to extend to infinity in the exterior domain.^i^ When the eigenmode (E, uj) is computed, the respective Q-factor 
is deduced from the real and imaginary parts of the complex eigenfrequency, Q = 5R(a;) / (— 23(w)), the resonance 
wavelength Aq is given by Aq — 27rco/w, with the speed of light cq. 

Figure [S] shows the convergence of the computed Q-factors and resonance wavelengths for the cases inves- 
tigated also in Section [5711 {N = 2/3/4). Here we have increased the numerical resolution by increasing the 
polynomial degree p of the chosen finite element ansatz-functions. As can be seen from this figure, both, Q- 
factor and resonance wavelength computed with the eigenmode solver converge to the values obtained from the 
light scattering simulations (c./. Section [O]) . Figure H] shows the computational costs corresponding to these 
simulations. Please note that simulating a resonance directly requires a single computation only while deducing 
the resonance properties from a transmission or reflection simulation (as in Section l3.ip always requires several 
computations at various wavelengths. Especially for high-Q resonances, where the choice of such wavelengths 
for a transmission scan is not clear a-priori, direct computation of resonances greatly simplifies the simulation 
task and reduces computational effort. 

3.4 Simulation of high-Q resonators 

In a 3D photonic crystal setup one expects that the cavity Q-factor will increase exponentially with the width of 
the surrounding photonic crystal when a cavity is operated inside a bandgap. However, in the given ID photonic 
crystal setup, the situation is different because leakage to the substrate and to surrounding air regions also have 
to be taken into account. In a further set of numerical experiments we have applied the eigenmode FEM solver 
in order to adress the question which Q-factors can be reached by simply adding more and more pores to the 
ID PhC microcavity. Figure [7] (left) shows that the cavity Q-factor increases exponentially with the size of the 
photonic crystal region for numbers of pores in the periodic sections of the design, N, between 1 and about 6. 
For higher N, a plateau of Q ^ 100,000 is reached. 

By integrating over the flux of energy density over the computational domain boundaries in the different 
material regions we have further recorded the losses to the substrate, to the surrounding air region, and to the 
waveguide region. Figure [7] (right) shows the corresponding results. In accordance with Fig. [7] (left), for TV 
below ~ 6, the losses to substrate and surrounding air region do not play a significant role (note the logarithmic 
scale of the plot). As in this regime the resonator looses its energy mainly through the waveguide, the photonic 




Figure 7. Left: Q-factor of PhC microcavities in dependence on number of air pores in the periodic sections of the design, 
A^. Right: Normahzed energy flux to waveguide, substrate, air regions in dependence on number of air pores in the 
periodic sections of the design, A'^. 



crystal perforation of the waveguide yields an exponential decay of the field and exponentially lower losses with 
increasing the width of the waveguide. For iV > 6 light scattering to the substrate and to the surrounding air 
regions are the dominant loss channels, and further improvement of the photonic crystal in the waveguide region 
does therefore not lead to an improvement of the quality factor beyond the values reached at iV ^ 6. 

4. CONCLUSION 

Finite-element solvers have been used to numerically investigate photonic crystal microcavities. The results have 
been validated by comparison to experimental and theoretical results from the literature. Further, numerical 
solutions from a scattering solver and numerical solutions from an eigenmode solver have been shown to yield an 
excellent agreement of physical results within the limits of numerical errors. The influence of nanometer-scale 
geometry variations (placement and sidewall angle) on the performance of one-dimensional PhC microcavities 
has been investigated. Losses to substrate and air regions have been identified as limiting factors to ultra-high 
Q-factor performance of the investigated specific photonic crystal microcavity design. 
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